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Probing dark radiation with inflationary gravitational waves 
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Recent cosmological observations indicate the existence of extra light species, i.e., dark radiation. 
In this paper we show that signatures of the dark radiation are imprinted in the spectrum of 
inflationary gravitational waves. If the dark radiation is produced by the decay of a massive particle, 
high frequency mode of the gravitational waves are suppressed. In addition, due to the effect of the 
anisotropic stress caused by the dark radiation, a dip in the gravitational wave spectrum may show 
up at the frequency which enters the horizon at the time of the dark radiation production. Once 
the gravitational wave spectrum is experimentally studied in detail, we can infer the information on 
how and when the dark radiation was produced in the Universe. 



I. INTRODUCTION 

Recently, there are increasing evidence of the extra 
non-interacting relativistic degrees of freedom, in addi- 
tion to the standard three (nearly) massless neutrino 
species. The abundance of relativistic component is 
parametrized by the effective number of neutrino species, 
Neg, as 
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where p rc \ is the total relativistic energy density and 
p y = (ir 2 /15)T^ denotes the photon energy density mea- 
sured after the e + e~ annihilation with T 7 representing 
the photon temperature. The standard model predicts 
JVeff = 3.04. 

The N e s can be constrained from various observations. 
First, increasing N e s leads to larger Hubble expansion 
rate at the big bang nucleosynthesis epoch, which in 
turn results in increase of the primordial helium abun- 
dance. Recent observations suggest -/V c fF = 3.681q 70 at 
2a level (See, however, also Ref. Q for discussion on 
the error estimation in the helium abundance.) 

The cosmic microwave background (CMB) anisotropy 
is also sensitive to N g. The information on iV e fj is im- 
printed in the CMB anisotropy in some ways. First, in- 
crease of N c g makes the early integrated Sachs- Wolfe 
effect more efficient, and the first peak of the CMB 
power spectrum is enhanced. Second, it tends to make 
the scale of sound horizon smaller at the recombination 
epoch, resulting in shift of the peak positions in the CMB 
power spectrum toward high multipole moment. Third, 
it erases the small scale power spectrum due to the ef- 
fect of free-streaming. The WMAP seven-year results 
combined with standard rulers give N e s = 4.32^q|| at 
la level [3J. Adding small scale CMB measurements im- 
proves the accuracy as N e g = 4.56±0.75 for ACT Q and 
N cS = 3.86 ± 0.42 for SPT g| both at la level. Ref. @ 
combined the WMAP, ACT and SPT datasets with stan- 



dard rulers and obtained N e s — 4.08 



-0.71 

'-0.68 



at 2cr level.* 1 



See also recent related studies [1-[TT|. 

To summarize, at the current situation, observations 
suggest AN e ff = N c g — 3 ~ 1 at nearly 2a level. Mo- 
tivated by these increasing evidence of the extra light 
species, which is often called "dark radiation", models 
to explain dark radiation were proposed [12-26]. Al- 
though there are many candidates, if the dark radiation 
has only extremely weak interaction with the standard 
model particles, it may be difficult to detect it experi- 
mentally. Thus it is important to study how to confirm 
and distinguish models of dark radiation by other obser- 
vations.* 2 

In this paper we consider a novel method to detect 
dark radiation through inflationary gravitational waves 
(GWs). It is known that relativistic free streaming fluid 
can contribute to the anisotropic stress, which potentially 
affects the propagation of GWs [28j . This effect was con- 
cretely studied for the free-streaming neutrinos. GWs 
entering the horizon after the neutrino freezeout dissi- 
pate their energies and, as a result, a modulation feature 
shows up in the GW spectrum [28M30j j . It was also stud- 
ied in the context of large lepton asymmetry [3l| . There- 
fore, it is expected that the dark radiation also induces 
similar effects on GWs. As opposed to the case of neu- 
trinos, we do not know when and how the dark radiation 
was generated in the Universe. Thus the position and 
strength of the modulation in the GWs, if detected, tells 
us exactly about the production mechanism of dark radi- 
ation. In particular, models of dark radiation produced 
by decay of non-relativistic fields [T^ . fl5l . [l6l . [T8l - [23( | shows 
characteristic features in the primordial GW spectrum. 
The feature consists of combination of the anisotropic 
stress effect and the modified background expansion his- 
tory. This is detectable in future space-based GW detec- 
tors such as DECIGO 32] and BBO [H, HJ]. We also 
show that the GW spectrum will be a powerful tool for 



Note that the statistical significance depends on the prior on 



the Hubble parameter 0. 
# 2 In Refs. [lSL HH , the possibility that the dark radiation has 
(non-Gaussian) isocurvature perturbations was considered. In 
Ref. [27j the effect of dark radiation on CMB B-mode spectrum 
is discussed. These may provide additional information on the 
nature of dark radiation. 
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confirming the dark radiation produced thermally and 
decoupled at some epoch in the early Universe. 

This paper is organized as follows. In Sec. HI] we re- 
view a model of dark radiation produced by decaying 
particles. In Sec. IIIII we calculate the evolution of grav- 
itational waves in the presence of anisotropic stress in- 
duced by dark radiation, and show that characteristic 
signatures appear in the spectrum. Sec. HVl is devoted to 
conclusions and discussion. 



II. DARK RADIATION PRODUCTION BY 
DECAYING PARTICLES 
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A. Background evolution 

We consider the case where the non-relativistic matter 
4> decays into X particle which plays the role of dark 
radiation. Thus X is assumed to be massless and has 
no interaction with other fields. To be more precise, X 
must be relativistic until the recombination epoch and its 
interaction must be so weak that remains to be decoupled 
from thermal bath after the production by <f> decay. The 
evolution equations of components are given by 



P0 + 3i?p0 = -r^p^, 

Prad + 4i?p lad = T^l 

px + 4:Hp x = T^Bxp 



Bx)p<p, 



(2) 
(3) 
(4) 



where the dot represents time derivative, and the Fried- 
mann equation, 



H 2 



Ptot _ P<f, + Prad + PX 



3M 2 



3MJ, 



(5) 



where p$, p ra d and px are energy densities of (f>, visible 
radiation and dark radiation, respectively, Mp is the re- 
duced Planck scale, is the decay rate of <f>, and Bx 
denotes its branching fraction into X. 

The extra effective number of neutrino species is given 

by 



AN, 



off 
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10.75 " 
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(6) 



H«r* 



where g* s (T^) denotes the relativistic degrees of freedom 
at T = where the <j> decays, and px and p ra d are 
evaluated well after the <fi decay. In our numerical study, 
we take the standard-model value of g* s {T^,) = 106.75. 

In order to obtain AN e s — 1, the energy density of <f> 
should nearly dominate the Universe at the decay. There- 
fore, the expansion rate of the Universe around the </> de- 
cay epoch is modified. Fig. Q] shows the product tH as a 
function of cosmic time t normalized by idee, defined by 



tde 



1 



(7) 



Here we have fixed initial conditions of p$ and p ra d so 
that AN e s = 1 is realized. Solid (red), long-dashed 



FIG. 1: Evolution of the product tH as a function of cosmic 
time (normalized by t^cc) for Bx = 0.26 (red solid), 0.5 (green 
dashed), 0.7 (blue dotted) and 1.0 (magenta dot-dashed) for 
explaining AN c g = 1. 



(green), short-dashed (blue) and dotted (magenta) lines 
correspond to Bx = 0.26, 0.5, 0.7 and 1.0, respectively. 
It is seen that c/> has significant energy fraction around 
its decay and the expansion rate is modified from the 
radiation-dominated one (tH = 1/2). In the limit of cj) 
domination, we need Bx — 0.26 for A7V c ff = 1. Thus 
Bx > 0.26 is required in order to realize AN c g = 1. 
Since the background expansion rate is imprinted in the 
GW spectrum [35l - l50j . a particular shape in the GW 
spectrum is expected if dark radiation is produced by 
decaying matter, as we sill see. 



B. Model 

As one of the motivated models of <j> and X, we con- 
sider the saxion and axion in a supersymmetric axion 
model jjyl. This possibility was studied in Refs. [H, HE 
EE HE Ha] in the context of dark radiation. 

The axion is a pseudo Nambu-Goldstone boson as- 
sociated with the spontaneous breakdown of the global 
U(1)pq symmetry [52[. It solves the strong CP problem 
in the quantum chromodynamics. The axion has interac- 
tions suppressed by the U(1)pq breaking scale, f a . The 
value of f a is phenomenologically constrained as 10 9 GeV 
< fa < 10 12 GeV, and the axion mass is - 10~ 2 -10 -5 eV 
for this range of f a 53] . Thus the axion is a good candi- 
date of dark radiation. 

In a supersymmetric extension of the axion model, 
there appears a scalar partner of the axion, called saxion, 
which is massless in supersymmetric limit but obtains a 
mass from supersymmetry breaking effects. Writing the 
saxion mass as m^, the saxion decay rate into the axion 



3 



pair is typically [54 



1 rn% 
64tt p a 



(8) 



Assuming that the saxion decays in the radiation dom- 
inated era, the temperature at the saxion decay is esti- 
mated to be 



T^3x 10 b GeV 



( m <t> y 

V10 3 TcV/ 



3/2 /10 10 GeV 



fa 



(9) 



The saxion with mass of O(10 3 ) TeV is plausible by tak- 
ing account of the preference for high-supersymmetry 
breaking scale [55| in light of the recent discovery of the 
Higgs boson mass of 125 GeV (56[. The saxion often dom- 
inantly decays into the axion pair (Bx — 1)- The pro- 
duced axions are never thermalized below the tempera- 
ture - 10 7 GeV for f a > 10 10 GeV g7|. The abundance 
of relativistic axion after the (f> decay is then estimated 
to be 



Px 

Prad 



B 



X 



Ptot 



6 Td, \M P 



(10) 

where Tr is the reheating temperature after inflation and 
is the saxion initial amplitude. Therefore, for appro- 
priate choices of Tr and the axion abundance pro- 
duced by the saxion decay can account for the dark ra- 
diation : AN c g ~ 1. 



III. SPECTRUM OF GRAVITATIONAL WAVE 
BACKGROUND WITH DARK RADIATION 



A. Evolution equations 



where denotes the polarization tensor. As shown in 
Appendix, (t, k) satisfies the following equation 



+ 3Hh^ + !L h W 



2AH 2 



1 



a 4 (t)p tot (t) 



h w (t',-k.)dt', 
(13) 



where 



it-*- 



sin(u) 3cos(u) 3sin(u) 

— t 1 + — r 2 ' 14 

-iiO rti^i n,D 



with j2 being the second-order spherical Bessel function. 
Contrary to the case of neutrinos studied in Refs. [H,[13> 
Px(t') is inside the time integral since px does not scale 
as a while X is produced by the </> decay. In terms of 
U and u! defined as 



dt' 

*' dt" 



(15) 
(16) 



Eq. (fT3"|) becomes 



= - 24H , 



1 



u a A (u)p tot (u) 
x / a i {u')px{u')K(u — u') h' m (u',k)du', 



(17) 



where the prime denotes the derivative with respect to u 
and H u = a' /a. We have solved Eq. (fl7|) together with 
the background evolution @ - (H]) to derive the present 
GW spectrum. 



Now let us study the evolution of primordial GWs un- 
der the presence of dark radiation. The GW corresponds 
to the tensor perturbation of the metric. We define the 
line element as 



ds 1 



-dt 2 + a 2 (t)(<% + h l3 )dx l dx\ (11) 



where hij is the transverse and traceless part of the met- 
ric perturbation, and the Fourier amplitude of hij as 



f d 3 k 

hij(t,x) = J — fcyftky 1 " 



E 

A=+,x 



d 3 k 

(27r)» 



/i (A) (i,k) 



,(A) ikx 



(12) 



B. Overall normalization 

Before showing the detailed results, we here comment 
on the normalization of the present GW energy den- 
sity. During inflation, quantum fluctuations of the tensor 
perturbation is continuously generated which turn into 
stochastic GW background in the present Universe af- 
ter the horizon- in [58|. It predicts nearly scale invariant 
GW spectrum for the GW modes entering in the horizon 
in the radiation-dominated era [59J-166J . The GW energy 
density per log frequency at the horizon crossing k = aH, 
normalized by the critical energy density, is given by (oTj 



^Gw(fc = o,H) = 



A?0) 2.43 x 10- 9 r 



24 



24 



TT ,(18) 



4 



where r denotes the tensor-to-scalar ratio, n* is the tensor 
spectral index, fc — 0.002 Mpc -1 is the pivot scale and 



A*(*) - Wp 



H 



inf 



2tt 



fc , 



(19) 



with -ffi n f being the Hubble scale during inflation and we 
have assumed the WMAP normalization on the curva- 
ture perturbation on large scale Q. In this subsection, 
we consider the modes which enter the horizon in the 
radiation-dominated era, since we are interested in the 
high-frequency GWs which may be observed by space- 
based GW detectors. 

First, in the standard model without dark radiation, 
the present spectrum of GW is given by 

ng^(*) = 7 ^nHf x n GW (k = aH), (20) 



where Q 



(SM) 
rad 



4.2 x 10 h with h parameterizing the 



present Hubble parameter as Hq = 100/ikm/s/Mpc and 

4/3 
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(SM) 



g*(T in (k)) 

ASM) 

y*o 



ASM) 
1*s0 



g*s(T in (k)) 



(21) 



where g*Q M ^ = 3.36 and = 3.91, and T- m (k) denotes 

the temperature at which the mode k enters the hori- 
zon. We have 7< SM ) ~ 0.39 for g«(T in (fc)) = 106.75. The 
present GW spectrum per log frequency is then given by 



n 



(SM) 
GW 



(k) ~ 3.3 x 10~ 16 



(o.l) ( fc 



106.75 



g*(T in (k)) 



1/3 



.(22) 



In the presence of dark radiation, the overall normal- 
ization of the GW spectrum is modified due to the change 
of expansion rate. Neglecting the effect of anisotropic 
stress, we find 



^Gw(fc) = 7^rad x Q GW (k = aH), 
where fi rad = O,^' x (#*o/s*o ) witn 



.9*0 — 2 



1 + iV, 



err: 



8 V 11 



4/3' 



(23) 



(24) 



We find 5*0 — 3.82 for N e g = 4. The factor 7 is given by 

AN cS 

- = — -T/3 > (25) 



1 



^ ^ 43 \ 10.75 J 

/7 (sm) + r 



AJV n1 



where we have used the relation ([6]) . Therefore, the over- 
all enhancement factor for the GW spectrum is given by 



C x = 



7 



5*o 



~(SM) (SM) ' 

' 9*0 



(26) 



o 




FIG. 2: Relative normalization of the GW spectrum for 
few <C k -C fcdec as a function of AiV e ff. Ci (green dashed) : 
the enhancement factor due to the modified background evo- 
lution. C2 (blue dotted) : the suppression factor due to the 
anisotropic stress caused by dark radiation. Their product 
Ci x C2 (red solid) determines the resulting overall normal- 
ization for the modes &ew "C k <C fcdec- 



which is 1.35 for AN c g — 1. The first factor comes from 
the modified expansion rate between the horizon-in and 
matter-radiation equality, and the second factor comes 
from the change of the epoch of matter-radiation equal- 
ity. Thus, without the effect of anisotropic stress, the 
GW amplitudes at high frequencies inferred from the 
measured tensor-to-scalar ratio at the CMB scales are 
enhanced in the presence of dark radiation. 

Such an enhancement is compensated by the dissipa- 
tion of GWs caused by the anisotropic stress of dark ra- 
diation. The suppression factor due to the anisotropic 
stress, which we express here by C2, was analytically de- 
rived in Ref. [2!| [38| as a function of energy fraction of 
relativistic free-streaming particles. The relative normal- 
ization for the GW spectrum is then given by the product 
of them, 



fi G w(fc) 



(SM) 
GW 



(A;) 



Ct x C 2 , 



(27) 



for k &dcc7 where fcdec denotes the comoving Hubble 



a(tdoc)H(tdec)- 



(28) 



Fig. [2] shows C\ , C2 and their product as functions of 
AiV e )f for A;ew "C k <C fcdec where A;ew denotes the co- 
moving Hubble scale around the electroweak phase tran- 
sition. It is seen that there is a cancellation between C\ 
and C2, and the result is close to one for AN c s = O(l). 
Although Eq. (|2"7| gives normalization of the GW spec- 
trum for /cew "C k -C fcdoc, the precise shape of the 
GW spectrum around k <~ fcdec needs to be investigated 
numerically. Detailed results are shown in the next sub- 
section. 
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FIG. 3: f2Gw(&)/^GW W as a function of fc (normalized by 
fcdec) for = 0.26. The green (dotted) line is the full result, 
taking account of the effect of anisotropic stress. For compar- 
ison, in the solid line (red), we also plot flGw(k) /Oq^ 1 (k) 
without the effect of anisotropic stress. 



C. Results 

In Figs. [3] - El we plot the GW spectrum normalized 
by Oq-vv CO predicted in the present scenario, varying 
Bx from 0.26 to 1.0. The horizontal axis is normalized 
by k dcc . For comparison, we have also plotted the GW 
spectrum without the effect of anisotropic stress. As one 
can see, the spectrum of the GWs has a characteristic 
change at k ~ fcdec if the dark radiation (with AN e g ~ 
1) is produced by the decay of massive particle. Thus, 
once the GW spectrum is precisely measured, we have a 
chance to extract the information on the mechanism of 
dark-radiation production. 

There are several effects on the GW spectrum in the 
presence of dark radiation. First, since 4> (nearly) dom- 
inates the Universe at the decay in order to realize 
AN e ff ~ 1, f^GW decreases at k > fcdec- This is due to the 
change of equation of state of the Universe. As a result, 
as one can see, JJqw is suppressed for high frequency 
modes which enter the horizon before the (^-domination. 

In addition, most importantly, the effect of anisotropic 
stress caused by dark radiation dissipates the GW en- 
ergy density of the mode with k < fcdec ■ Consequently, 
together with the effect of the change of equation of state, 
a dip in the spectrum may appear at k ~ fcdec- In partic- 
ular, the dip becomes more apparent when Bx is close 
to 1 . Such a dip provides a smoking-gun signature of the 
dark-radiation production by the decay of massive par- 
ticles. If 4> and X are completely sequestered from the 
standard-model sector, for example, Bx = 1 may be re- 
alized. Then, such a model provides a striking signature 
in the GW spectrum. 

Note that, in the low frequency limit k <C fcdec j we have 



FIG. 4: Same as Fig. ED but for B x = 0.5. 
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FIG. 5: Same as Fig. E3 but for Bx = 0.7. 



numerically confirmed the suppression factor Ci caused 

by dark radiation. As a result, ficw/^Gw'' at fc <C fcdec 
is close to one as shown in Fig. [2] 



IV. CONCLUSIONS AND DISCUSSION 

In this paper we have studied the spectrum of infla- 
tionary GW background in the presence of dark radia- 
tion, motivated by recent observational preferences for 
AiVeff ~ 1. We have assumed that the dark radiation is 
non-thermally produced by decay of massive particles </>. 
There are several effects on the GW spectrum. First, the 
equation of state of the Universe is modified due to the <fi 
energy density and it changes the shape of the GW spec- 
trum. Second, the anisotropic stress carried by dark ra- 
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FIG. 6: Same as Fig. H but for B x = 1.0. 



the dip in the GW spectrum cannot be seen in the GW 
detectors. Instead, overall normalization of the GW spec- 
trum at the observable frequency range, inferred from the 
measured tensor-to-scalar ratio, is enhanced by the fac- 
tor C\ r~j 1.3. (At this epoch, dark radiation took part 
in thermal bath and there is no anisotropic stress damp- 
ing on GW amplitudes with corresponding modes.) This 
provides another indirect evidence of dark radiation. 
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diation dissipates the GW amplitude for modes entering 
the horizon around and after <f> decay. Numerical results 
show that there may appear a characteristic dip around 
k ~ fcdec, which is a smoking-gun signature of dark radi- 
ation. It not only provides an evidence of dark radiation, 
but also sheds light on its production mechanism. 

Some notes are in order. We have assumed that the 
dark radiation anisotropic stress is induced only by the 
primordial GWs. This is not in general true in the sec- 
ond order perturbation theory. Free-streaming particles 
(as well as other fluids) contribute to GWs at the second 
order in the scalar perturbation even if there is no pri- 
mordial tensor perturbation. However, this contribution 
is negligible for r > 10" 6 [Hi!]. 

So far, we have considered dark radiation produced by 
the decay of <j>. However, it is possible that the dark 
radiation was once in thermal equilibrium and decoupled 
from thermal bath at the temperature Tdcc- In this case, 
the extra effective number of neutrino species is given by 



AN, 



off 



-eN x 



10.75 

Q*s (^dec) 



4/3 



where 



for a real scalar, 
for a chiral fermion, 



(29) 



(30) 



and Nx counts the number of X species. If the de- 
coupling temperature is higher than the weak scale, 
we need Nx > 20 for explaining AiV c ff ~ 1. The 
modulation in the GW spectrum, similar to the effect 
caused by of neutrinos apparent at the GW frequency of 
10~ 10 Hz 0, appears at the frequency inside the 
range of DECIGO/BBO sensitivities for T dcc - 10 7 - 
10 9 GeV. If the decoupling temperature is 0(1) MeV, 
Nx ~ 1 is sufficient in order to obtain AN e g ~ 1 but 



Appendix A: Equation of motion of gravitational 
waves with dark radiation 



In this Appendix we derive Eq. (|17[) . the equation 
of motion of GWs with dark radiation. We follow 
Refs. [28|, |30( but the result is slightly different because 
X is continuously produced by the decay of <f> so that the 
number of X in the comoving volume is not constant. 

Throughout this appendix, we use the synchronous 
gauge and consider tensor perturbations defined in 

Eq. mj. 

The equation of motion for tensor perturbations in 
Fourier space is 



hij = 16ttG (~) n 



(Al) 



where the prime denotes the derivative with respect to 
u = k J dt/a and is defined by using the total energy 
momentum tensor as 



T, 



(tot) 



p 



Pg 



(tot) 



(A2) 
(A3) 



Our goal in this appendix is to express the RHS of 
Eq. (|A1[) in terms of metric perturbations. In what fol- 
lows, use the fact that only the collisionless particle (i.e., 
X) contributes to the anisotropic stress Ily. 

We first introduce the distribution function of the rela- 
tivistic components F' tot ' (t, x l ,pi), with which the total 
number of relativistic particles with particular momen- 
tum range contained in the volume element is given by 
F( to Vdx 1 dx 2 dx 3 dpidp2dp3. (Here and hereafter, x % is the 
comoving coordinate, while pi is the comoving momen- 
tum.) Note that F' tot ) is a scalar under general coor- 
dinate transformations which preserve the synchronous 
gauge. The distribution function can be decomposed as 

F( tot ' (t,x i ,p i ) = F( x \t, x l ,pi) + (t,x l , Pt ), ( A4) 
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where F( x ) and F( rad ' are distribution functions of the 
dark radiation X and that of ordinary radiation (like 
photon, gluon, and so on) with very short free-streaming 
length, respectively. Hereafter, we omit the superscript 
X for the distribution function of X for notational sim- 
plicity : F{t,x\ Pi ) = F^ x \t,x\ Pi ). 

We start with the effect of dark radiation on the 
anisotropic stress. The distribution function of X obeys 
the collisionless Boltzmann equation with source from 
non-relativistic decaying particle 4>: 



dF B x / 



4tt(p°) 



2 ) 



(A5) 



where p° is the energy of X , and we assume that <fi decays 
into two Xs. Also note that p° and p l should be regarded 
as functions of pi through g^vP^p" — — (p ) 2 + a 2 (<% + 



hij)p l p J = and p l 



of Eq. flA5]) is 



9 l3 Pj 



' (Sij — hij )pj . The LHS 



dF dF dx l dF dpi dF 

= 1 7 + — 

dt dt dt dx l dt dpi 

dF p i dF t 1 dF 
p° dx % 



(A6) 



dt 



d Pk 



where we used 



dx l p 1 

~dt = p°' 

dpi 1 pPp k 

~dt = 2 93k ' l ~p^- 



(A7) 
(A8) 



Eq. (|A8|) is obtained from the geodesic equation. 

Next we decompose F into the unperturbed part 
F(t,p), where p = ^Jpipl should not be confused with 
the pressure, and the perturbed part SF. We further 
decompose SF into two terms SF\ and 6F 2 for later con- 
venience: 



5F 1 {t,x i ,p i ) = FfatftpiptfP/at-Ffap) 

SF 2 (t,x i ,p i ) = F-F-SFl 
We get from Eq. (|A5I) the zeroth-order equation 
dF 



(A9) 
(A10) 



It = a^y t ^ 5 ( p0 ~ ' (A11) where 



and the first-order one 

d{5F 1 + 5F 2 ) p 1 d(SF 1 + 5F 2 ) 



dt 



1 



2 -=o 



p° dx % 



p° dpi 



dpdt 



Sp°. 



(A12) 



In Fig. [71 we show 4irp 3 F as a function of p/pdeo where 
Pdcc = a{tdec)fn^/2 is the comoving momentum of X 
produced ait = idec- We can see that the energy fraction 



rd 
o 

co >, 
CO 



100 



10 



I I I I 1 — 

t=io 4 w 



- 1 1 — i — m 



- 1 1 1 - 




P/Pdec 



FIG. 7: Aitp 3 F as functions of p (normalized by pdec = 
a(tdec)m t / > /2). Each line corresponds to t > t^cc (red solid), 
green long-dashed), 10 2 idec (blue short- 
dotted) and tdec (magenta dotted). 



4irp 3 F is mostly carried by X produced at t ~ £dcc- Then 
we use the following equations: 



SF X = --hij 



5 P ° = 



l u PiPjdF 
p dp ' 
1 hijpiPj 



2a p 



(A13) 
(A14) 
(A15) 



Using p t =Pi/ p and substituting Eq. (|A13|) - Eq. (IA15|) 
into Eq. (|A12I) . we get 

dSF 2 p % dSF 2 _ 1 dhij dF 



+ a^ = 2^T^^- (M6) 

In terms of conformal time t — J dt/a, this equation is 
expressed as 



dSF 2 



P 



dr dx % 
In Fourier space, 



dSF 2 _ 1 dhij dF 
2 dr dp 



PPiPj- 



dSF 2 
dr 



ikfiSF'^ 



1 dhij dF 



2 dr dp 
d 3 k 



PPiP v 



(A17) 



(A18) 



SF 2 (r,x\ Pi ) = J JL^6F 2 (T,ki,pi)e ikixi ,(Al9) 

i r i\ f d 3 k 



rhij(T, h)e 



(A20) 
(A21) 



(2tt) 3 
P = hpi. 

We can use line-of-sight integral to get the solution of 
Eq. (|AT8|t : 

tf2= l d 4¥ (T ' ) S (T ' )P « e " ,l '' lM) ' (A22) 
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where we have used SF2(t — 0) = because there is no 
X in the beginning. 

We take the first-order perturbation of the energy- 
momentum tensor of X: 



T { * ] = 1 



(A23) 



(A24) 



Note that energy momentum tensor defined above trans- 
forms as a tensor under general coordinate transforma- 
tions s ince / d 3 p/p° oc d A p5{g^ u p^p v ). Using Eq. (|A13|) 
- Eq. flMgJ) , we get 



<^ ( f 5 = ^ I d 3 p SF 2 app i p J 



1 , 2 . „ „ „ 1 , . . A „ - 

--ahkip PiPjPkPi-^ + ^ah k ippiPjp k piF 



1 



d 3 pSF 2 appiPj 



dpp 



1 , ,9F 1 , 
-ah u p — + -ahkipF 

2 op 2 



47T 



If 1 

= -2 / d 3 p6F 2 pp l p j + -a 2 h ljPx . (A25) 



Here, we used 

dtt p piPjpkPie~ tpii:iU 

47T 

, J2(«) 



ji(u) ~kikjk k ki - J ' 3 ^ (hkjSki + 5 perms) 



and 



-(SijSki + 2 perms) 



dp4np F — a px, 



(A26) 



(A27) 



where is the energy density of X and j n is the n-th 
spherical Bessel function. 

Next, we consider the effect of _F( rad ) ; for which 
5^2 = because the free-streaming length is very 



short. Then, we obtain 



rrTi(rad) J- 2 u 



(A28) 



We also note that perturbation in the energy momentum 
tensor of <j> vanishes since it behaves as non-relativistic 
matter : 

ST^ = 0. (A29) 
Taking the first-order perturbation of Eq. (|A3|) . we ob- 
tain 



5T {t0t) = SP . g.. + p . 5g .. + a 2 U .. 

= ^a 2 hij(px + Prad) + a 2 nij, (A30) 



wher e we used Eq. (jA"T3]) - Eq. (|AT5l) . Eq. (|A"22l . 
Eq. (|A26)l . h« = and SP = 0. The last condi- 
tion comes from the fact that tensor perturbations can- 
not produce perturbations in scalar variables. Using 
Eq. (|A"25l> . Eq. (|A"28l) . Eq. (|Al9|l , and Eq. (|A"30)l . we ob- 
tain 

a 2 n„ = i /" d 3 P SF 2 pp lPj . (A31) 



Substituting Eq. (|A"31~]) into the RHS of Eq. (fATj) . we 
obtain 



/a\ 2 1 /" 
= 16ttG (-J -j / d 3 pppipj 

r , ,idh kl( ,.dF , 



k 2 a 2 



I 



Or 



dClp— a pxPiPjPkPie 



— ikfl(T-T') 



= -24H„ — — / duapx-^-iu) 



«„4 



a 4 /9tot(w) Jo 



9m (m — ?i') 2 

(A32) 



where we used partial integration, Eq. (|A26j) and Fried- 
mann equation i/ 2 = 87rGptota 2 /3fc 2 . After decomposing 
using Eq. (fl"2j) . we finally obtain Eq. (fTTl) . 
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